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We present in this paper a compact model of bipolar transistors, suitable 
for network analysis computer programs. Through the use of a new charge 
control relation linking junction voltages, collector current, and base charge, 
the model includes high injection effects. The performance substantially 
exceeds that of existing models of comparable complexity. For low bias 
and with some additional idealization, the model reduces to the conventional 
Ebers-Moll model. 

I. INTRODUCTION 

Since its formulation in 1954, the Ebers-Moll model 1 has been the 
major large-signal model for bipolar transistors. It is based directly 
on device physics and covers all operating regimes, that is, active, 
saturated, and cut-off operation. But various approximations limit the 
accuracy of the model. In the original paper the frequency dependence 
was described in terms of frequency-dependent current generators. 
For time-domain analysis, this required use of Laplace transforms. 

In 1957 Beaufoy and Sparkes 2 analyzed the bipolar transistor from a 
charge control point of view. The charge control model 2 or the equiva- 
lent 3,4 charge control form of the Ebers-Moll model, is directly useful 
for transient analysis. It is this model, in various forms, that is presently 
used in all the major general purpose network analysis programs. More 
elaborate models 5-7 are useful for detailed device studies, but they have 
not found widespread use in analysis programs due to their complexity. 

As device technology evolved over the years making possible devices 
of reproducible characteristics, and as deeper understanding was gained 
in device-theoretical studies, many new effects were identified that are 
not represented by the basic Ebers-Moll model. Among these are a 
finite, collector-current-dependent output conductance due to basewidth 
modulation 8 (Early effect), space-charge-layer generation and recombi- 
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nation 9 (Sah-Noyce-Shockley effect), conductivity modulation in the 
base 10 (Webster effect) and in the collector 11 (Kirk effect), and emitter 
crowding. 13 ' 18 

Some of these effects have been included into the Ebers-Moll model 
in various ways. The common approach has been to specify some param- 
eters of the model as functions of bias and to describe this bias de- 
pendence in tabular form or through parametric expressions. Thus, 
in the NET-1 program, for example, 14 the common-emitter current 
gain is given by a series expansion in the emitter-base voltage. In the 
CIRCUS program, 15 forward and reverse current gains and forward 
and reverse transit times are specified as functions of collector current 
in tabular form. Such "curve-fitting" modeling tends to require large 
numbers of parameters or table entries for an accurate description. 
Also, frequently the parameters are not easily interpretable in terms 
of the device structure and thus can be obtained only a posteriori, from 
detailed measurements, but cannot be predicted. 

The present paper makes use of a new, general charge-control re- 
lation 18 which links junction voltages, collector-current, and base 
charge. This new charge-control relation, used in conjunction with 
conventional charge-control theory, allows many of the effects not 
contained in the basic Ebers-Moll model to be incorporated in an 
integral, physical way, in compact form and with good parameter 
economy. 

In the conventional treatment of the transport equation, as used in 
the derivation of the Ebers-Moll model, attention is focused on re- 
combination in the base (which in modern transistors is generally not 
the dominant source of base current). Through the use of simplified 
"boundary conditions" relating carrier concentrations at the edges of 
the base to the junction voltages, and through neglect of the variation 
of the electric field in the base with bias, the voltage dependence of 
emitter and collector currents contains errors, except near zero bias. 

The new charge control relation arises from the treatment of the 
transport equation for the carriers that pass between emitter and 
collector. Use is made of the fact that recombination has only a very 
small effect on the junction-voltage dependence of the current passing 
from emitter to collector (later called the dominant current component). 
Hence for this dependence, but of course not for the base current, 
recombination is neglected. A direct, closed-form solution of the trans- 
port equation from inside the emitter to inside the collector is then 
possible. An important feature of the result, quoted in equation (10), 
is a denominator containing the base charge, that is the charge of all 
those carriers that communicate with the base terminal. 
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With the concentration in the base of the carriers carrying the domi- 
nant current component known, recombination in the base could be 
computed, and the small effect that this recombination has on the 
dominant current component could be applied as a correction. 17 Actu- 
ally, however, for typical transistors this correction is negligible. The 
base current is modeled by empirical relations since detailed recombi- 
nation parameters are not available. 

The general features of the new model are presented in Section II. 
Details of a particular implementation using 21 model parameters are 
given in Section III.* Table I contains a list of the model parameters. 
Appendix A describes the treatment of base push-out. Discussion of 
some aspects of the model, and a set of computed characteristics of a 
sample transistor are given in Section IV. A compact summary of the 
model is given in Appendix B. A discussion of the significance of the 
model parameters and of default values is contained in Appendix C. 
This paper considers mainly the dc and low-frequency aspects of the 
model. Use of the model for ac or transient analysis will be presented 
elsewhere. 

II. GENERAL CONSIDERATIONS 

As a point of departure we will consider the Ebers-Moll equations. 1 
They may be written in the form 

?.\ = [T] /exp(qV tb /kT) - l\ (]) 



Jj \exp(qV cb /kT) - 1, 

where T is a symmetric matrix of coefficients that are constant, that is, 
bias independent. At the time when the Ebers-Moll model was de- 
veloped, attainable base widths were large by today's standards, and 
in order that useful current gains could be obtained lifetime in the 
device had to be long. Reverse saturation currents were used as indi- 
cators of lifetime. These circumstances are reflected in the notation 
used in Ref. 1 for the elements of T: 



T = 



I.o aj. 



1 — a„ai 1 — a n ai 

aj. I„ 



(2) 



1 — a n oti 1 — <x n ai 

where I.„ and I co are emitter and collector reverse saturation currents 

* A simplified implementation of model, and a comparison with measurements, 
have been presented in the Digest of Technical Papers, 1970 IEEE International 
Solid State Circuits Conference, Paper 7.1. 
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Table I — Model Parameters 

Group 1 : Knee parameters and transit times 
Ik Knee current* 

Vk Abs. value of knee voltage, in unit of kT/q 
Tf Forward tau (forward transit time) 
r t Tau ratio (ratio of reverse to forward transit time) 

Group 2: Base Current 

ii Ideal base current coefficient 

ij Nonideal base current coefficient 

n, Forward base current emission coefficient 

t*a Reverse base current coefficient 

■tic Reverse base current emission coefficient 

Group 3: Emitter Capacitance 

rv„ c Abs. value of emitter offset voltage, in units of kT/q 

\m e Emitter grading coefficient 
PAciei Emitter zero bias capacitance coefficient 

[a e 2 Emitter peak capacitance coefficient 

Group J f : Collector Capacitance 

rv oc Abs. value of collector offset voltage, in units of kT/q 
\m c Collector grading coefficient 
P c~\o, c \ Collector zero bias capacitance coefficient 
[flei Collector peak capacitance coefficient 

Group 5: Base Push-out 

v rp Abs. value of base push-out reference voltage, in units of kT/q 

r w Effective base width ratio 

r p Base push-out transition coefficient 

rip Base push-out exponent 

Auxiliary Quantities 

e k = exp ( — v k ) 
e ke = exp ( —Vk/Ue) 
eke = exp ( —Vk/ric) 

* Quantity is negative for npn transistor. 



and a n and a, are forward and reverse common base current gains. A 
simpler and more appropriate but fully equivalent notation is 

T = /(i + W,)i. -I. \ (3) 

\ - /. (1 + l/ff r )lj 

Here /3/ and /3 r are forward and reverse common emitter current gain 
and /, is the "intercept" current, that is, the current obtained when 
on a semilog plot of I c vs. V tb the collector current is extrapolated to 
V eb = 0. The notation of equation (3) is considered more appropriate 
since the intercept current I, is nearly independent of current gains. 
The matrix elements T l2 and T 21 in equation (2) show an apparent 
dependence on lifetime through the forward and reverse current gains. 
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Actually, this dependence is nearly cancelled by that of /,„ and I co . 
To a very good approximation the intercept current J, depends only on 
the total number of impurities in the base. 18 

Equation (1) with matrix T given by equation (3) suggests the 
following interpretation: The emitter and collector current have a 
common, dominant component 

I cc = -J.[exp (qVJkT) - cxp (qV ab /kT)]. (4) 

In addition the emitter and collector currents have each a separate 
component I be and I bc , which is proportional to the reciprocal forward 
and reverse common emitter current gain: 

I h . = (/.//3,)[exp (qVJkT) - 1], (5) 

he = (L/0r)[exp (qVJkT) - 1]. (0) 

The terminal currents are then given by 

/„ = -I ee + I be , (7) 

L = I cc + I be , (8) 

I b = -I be - I bc . (9) 

So far, we have made no changes in the content of the Ebers-Moll 
equations; we have only brought them into a form that will facilitate 
the development of the new model. We shall retain equation (7) through 
(9)* but we will replace equations (4) through (6) by relations giving 
an improved representation of the physical processes in the transistor. 

The separation of emitter and collector currents into the dominant 
I ee component and the base current components in equations (7) and 
(8) allows us to give different voltage dependence to the individual 
components. For example, at low injection levels, collector current 
and emitter-base voltage are related through the "ideal" diode law; 
that is, the collector current is proportional to exp (qV eb /nkT) where 
the "emission coefficient" n is very close to unity. The base current at 
low forward bias, on the other hand, is typically a "nonideal" current, 1 
that is it has an emission coefficient n with values typically between 
1.5 and 2. This nonideal current results from space charge recombina- 
tion or surface recombination, or the presence of both effects. At 
higher forward emitter-base voltages the base current is dominated 
by an ideal component. 10 

In principle, it is possible to compute the base current as a function 



These equations are not independent. Equation (7) will be dropped later. 
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of V. b (for given V eb ) for one-dimensional structures, 17,20 " 22 provided 
that the doping profile and the recombination parameters, for example 
the concentration as a function of distance of the important species of 
recombination centers, as well as their energy levels and capture cross 
sections, are known. However, in practice, the recombination properties 
are not known to the detail required for such calculations. Even the 
assumption of a constant concentration of one species of recombination 
centers is a gross oversimplification. In real transistors the lack of 
lattice perfection in heavily doped regions causes enhanced recombi- 
nation, and interfaces between substrates and epitaxial layers provide 
local regions of high recombination. 

Nevertheless, detailed studies have confirmed that the base current 
can be described by a sum of terms exponential in voltage with emission 
coefficients of the magnitude indicated above. Equations (5) and (6) 
are thus replaced by more general functions of V. b and V eb that are 
characterized by pre-exponential coefficients and emission coefficients 
which become model parameters and which represent the overall 
recombination properties and influence the dependence of forward 
and reverse current gain on bias. In the next section we will use a 
specific set of base current components. 

The new charge control relation 10 , derived from basic physical 
considerations, permits the dominant current components I ee to be 
obtained without use of low-injection approximations: 

T _ _ IQb aw (QV. JW) ~ expfoWttT) (1Q) 

Q b 

Here Q b is the "base charge", that is, the charge of all carriers of the 
type that communicate with the base terminal: electrons in a pnp 
transistor and holes in an npn transistor. The base charge is a function 
of bias; Q b0 is the zero-bias value. The bias dependence of Q b is the 
subject of the conventional charge control theory 12 which approximates 
the excess (or stored) base charge Q b — Q bo as the sum of capacitive 
contributions and products of collector current components and transit 
times. The essential feature of the model presented here is the simul- 
taneous use of conventional charge-control theory and relation (10). 
As shown before, 18 relation (10) is a good approximation for one- 
dimensional transistors for a wide range of bias conditions, including 
those leading to high-injection and base push-out effects. 

It is of interest to note that the Ebers-Moll equations embody 
superposition, that is that the collector current can be expressed as the 
sum of a function of the emitter voltage and a function of the collector 
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voltage. For real transistors, violations of the superposition principle 
are easily observed. Consider, for example, the Early effect 8 that is, 
the dependence of the low-frequency output conductance on bias. As 
shown schematically in Fig. 1, a region of bias exists in which the 
collector current varies approximately linearly with collector-emitter 
voltage for fixed base current, in such a way that the straight-line 
sections, when extrapolated, intersect (approximately) at a negative 
voltage which we shall call the "Early voltage", V A . For superposition 
to be valid, the lines would have to be nearly parallel to each other. 
The base charge Q b in the denominator of equation (10) through its 
dependence on collector voltage via the collector capacitance disables 
superposition and provides a realistic description of the output con- 
ductance. 

Another point of interest concerns high-injection effects in the base 
region. The "ideal" voltage dependence of I cc on V tb is caused primarily 
by the dependence of the minority carrier concentration in the base 
near the emitter as exp (qV. b /kT). This dependence holds, however, 
only as long as the minority carrier concentration is small compared 
to the doping concentration. If the minority carrier (the word minority 
starts to lose its literal meaning here) concentration is large compared 
to the doping concentration, then it varies as exp (qV. b /2kT), and so 
does, approximately, I ee except for additional complications due to 
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Fig. 1 — Definition of Early voltage. 
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base push-out. The intersection of the n = 1 and n = 2 asymptotes 
to I ee represents an important characteristic feature of the transistor. 
We shall call it the "knee point" and use its coordinates, V k and 7* as 
model parameters and as a basis for normalizations.* 

Consolidating the development up to this point, the model, exclusive 
of parasitic effects, is described by 

L = -LQ bo ex P faWftH - exp {qVJkT) + ^ (n) 

Qb 

/, = -I be - I bc . (12) 

As discussed above, the base current components I be and I be depend 
strongly on the recombination properties of the structure and are in 
practice not readily calculated from first principles. By contrast, the 
base charge as a function of bias depends primarily on the doping pro- 
file and is nearly independent of recombination properties. Hence, given 
the doping profile, Q b as a function of V eb and V cb can be computed by 
existing techniques. However, for such a calculation considerable com- 
puter resources (memory and time) are required. For network analysis 
programs it is preferred to approximate Q b by simple algebraic or 
algorithmic (implicit functions) representations which, depending on 
complexity, can give reasonable accuracy. One such representation is 
given in the next section. Special features are use of a modified rep- 
resentation of junction capacitance 23 which avoids the problem of an 
infinite capacitance when the junction voltage equals the built-in 
voltage, and use of a 4-parameter representation of base push-out. 

III. SPECIFIC MODEL IMPLEMENTATION 

In this section we present a specific implementation of the general 
model, equations (11) and (12). The bias dependence of base charge 
and base current will be modeled. The polarity assumed is that of a 
pnp transistor. 

The dominant current component I cc may be separated into an 
emitter and a collector component, or a forward and reverse component: 

T _ rn exp {gVJkT) - 1 

+ LQ ba exp (qV - /kT) - ' --/, + /,. (13) 



* The exact definition of these parameters is given in the following Section. 



CHARGE CONTROL MODEL 835 

We then model the excess base charge as consisting of emitter and col- 
lector capacitive contributions Q e and Q c , and of forward and reverse 
current-controlled contributions* 

Q b = Q ho + Q c + Qc - T f BI, - r,I r . (14) 

Here 77 and r r are forward and reverse transit times. The coefficient B 
has been included to model the increase of the transit time when base 
push-out occurs; it has a value of unity in the absence of base push-out. 
For modeling these parameters it is convenient to normalize all 
charges in equation (14) with respect to the zero bias charge Q bo and 
to denote the normalized charges by lowercase symbols. Also, we replace 
i/ and I r according to equation (13). Then 

1 + q, + q, + ~ -7T ! T / fi [ cx P (ffV".i/*D - 1] + r r [exp (qV eb /kT) - 1]} 

' ■ Qb ~ Vfto 

. 1 

ffi + - «■ ■ 

qb 

Qi 

Multiplication of equation (15) by q b removes q b from the denominator 
of the last term on the right side of equation (15) and gives rise to a 
quadratic equation in q b . Its solution gives q b explicitly in terms of the 
junction voltages, except for a possible g 6 -dependence of B: 



*-f + 



f) +* 



(16) 



The term q t represents the sum of the zero-bias charge and the charge 
associated with the junction capacitances; q 2 represents the excess base 
charge, or the current-dependent charge associated with diffusion capaci- 
tances. The latter charge contains a dependence on the base push-out 
effect through the parameter B which is treated in Appendix A. 

For high forward bias the charge q 2 is the dominant component of 
the base charge q h . Except for the base push-out term B, it is char- 
acterized by four parameters: /, , Q bo , r, , and r r . It will be convenient 
to normalize these parameters. For this we define the knee voltage V k 
as the emitter voltage for which q 2 equals unity (for zero collector voltage 
and with neglect of terms small compared to the exponential of the 
emitter voltage) : 



* The minus sign arises because the base charge contains electrons neutralizing 
the positive charges r f BI f and r T l T - Qb, Qbo, Qe, and Q c are all negative quantities 
for positive V eb and V C b- 
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7 A = ^ln-^- (17) 

The low-injection-extrapolated collector current for V eb = V k is 

h = — • (18) 

It will be convenient to normalize all quantities having dimensions of 
current with respect to 7* , and to express voltages by their difference 
from V k , in units of kT/q. We shall use lower-case symbols for normal- 
ized quantities. Thus, we define 

y * = W ' < 19 > 

Vt = 3 ^~ vk ' (20) 

v e = ^ - v k , (21) 

e* = exp (-v k ). (22) 

With these normalizations, equations (11) and (12) become 

I e = I k i c - J t (- e " ~ *" + 4 C ) , (23) 

7 6 = hi* = /.(-»,. - i be ). (24) 
For the base charge we obtain 

q s = B(e" - e k ) + (r r /r r )(c" - e*), (25) 

9, = fl x /2 + [( 9l /2) 2 + rf, (26) 

Q 6 = -/.Trf, . (27) 

With these normalizations, we can use the set It, , f* , r$ , r, = r f lr T as 
the four model parameters describing q 2 for the case B = I. These 
four parameters constitute Group 1 of the model parameters listed in 
Table 1. 

The charge contribution from the emitter and collector junction 
capacitances will be considered next. The conventional representation 
of junction capacitance is through an expression containing three 
parameters: 

c - w^SZzr (28) 
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The parameters are V bU iit-in (which for silicon is typically PH.7V), 
the grading coefficient m, and the constant in the numerator which can 
be related to the zero-bias capacitance. This expression causes diffi- 
culties when the junction voltage V approaches the built-in voltage and 
C goes to infinity. In a real transistor, of course, a finite amount of charge 
is stored for all bias conditions, and the derivatives of charge with 
respect to junction voltages are finite. As shown in detail in Ref. 23, 
equation (28) can be modified so as to be free of singularities by the 
introduction of a fourth parameter which relates to the forward-bias 
capacitance inferred from measurements of transit time vs. emitter 
current. For compactness of notation, we express the four parameters 
as elements pi , p 2 , P3 , and p 4 of a four-dimensional vector P (see 
Table I and Appendix C). If we now define a function 



Kv, P) = p 3 



■ + . ti 



fe-'M: 



(29) 



(i + Pi y 

LApi 
then the normalized emitter and collector charges are given by: 

"■ = i*w • *•) - (30) 

«■ - '(w 1 ■ p ') ' (31) 

where P. and P e are the four-parameter vectors describing the emitter 
and collector junctions, respectively. These two vectors constitute 
Groups 3 and 4 of the model parameters listed in Table I. As discussed 
in Appendix C, some of the parameter values can be estimated or ap- 
proximated in terms of other model parameters. In any case, these 
parameters are readily amenable to numerical evaluation from the 
device structure. 

As mentioned in the previous section, the recombination in transistors 
is best handled through a description of the base current as a sum of 
exponentials in the junction voltages. Pertinent model parameters are 
pre-exponential factors and emission coefficients. For typical transistors 
the forward base current is adequately described by two components, 
one ideal (n = 1) and the other nonideal (n = n e ). For the reverse base 
current a single nonideal {n = n e ) component is adequate. We define 

e k = exp (-v*), (32) 

e k . = exp (-v k /n.), (33) 
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e kc = exp (-v k /n c ). (34) 
Then we let 

i bc = i^e" — e k ) + ia(e" /n ' — e kg ), (35) 

i,,c = i*{e Vc/nc - e kc ). (36) 

The quantities U , i 2 , n e , i 3 , n e are the five parameters which char- 
acterize the recombination behavior of the transistor. They are listed 
as Group 2 in Table I. 

The last set (Group 5) of the model parameters listed in Table I 
describes the base push-out effect. For its description four parameters 
are required. The details of the base push-out modeling are given in 
Appendix A. 

We thus have a total of 21 parameters which describe the transistor 
model. For several of these parameters default values can be used as 
discussed in Appendix C. 

The following features are a consequence of the normalizations used: 

(i) I k is proportional to the emitter area. All other parameters are, 
to first order, independent of area. Area scaling (neglecting com- 
plications caused by lateral effects such as emitter crowding) is 
achieved simply by changing the value of I k . This feature is 
particularly convenient for integrated circuit work, where transis- 
tors on a given slice differ nominally only in their lateral dimensions. 
(ii) For pnp transistors all model parameters have positive numerical 
values. For npn transistors two changes are required: a) I k must 
be made a negative quantity; b) the Boltzmann voltage kT/q must 
be given a negative value (or the Boltzmann voltage is given the 
sign of I k for pnp's and npn's). When this is done, the polarity 
of terminal currents and voltages is in agreement with standard 
practice (currents positive if flowing into the device). 
(in) The offset voltages V ot and V oe used in modeling the capacitive 
charges (see Appendix C) are approximately proportional to the 
absolute temperature. Hence use of constant, that is, temperature 
independent, values for the normalized quantities v oe and v oe imple- 
ments automatically the temperature dependence of the offset 
voltages. The normalized knee voltage v k is not, to first order, 
temperature independent, but varies with temperature approxi- 
mately as 

«,CT) = v k (T ) + ^(^f- l) (37) 
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where T„ is a reference temperature (for example, room tempera- 
ture) and V g is the band-gap voltage (1.12 V for silicon). 

IV. DISCUSSIONS AND EXAMPLES 

This Section will illustrate the performance of the new model. The 
following features are of interest: 

(i) For low bias so that Q b is nearly equal to Q bo (or q b £rf 1) and with 

the choice n, — n e = 1, the model reduces to the Ebers-Moll 

model. 
(ii) Superposition which is operative in the Ebers-Moll equations, is 

disabled through the base-charge denominator Q b in equation (11). 

The dependence of Q b on collector voltage produces a finite output 

conductance, that is, the Early effect. 
(Hi) The rapid increase of Q b when base push-out occurs causes fall-off 

of current gain and frequency response. 
(iv) If wc define an effective emission coefficient n by 



1 = kT/q dl r 
n " I, dV„ 



(38) 



Fct-oon.t 

then it is seen that for the model, n varies from approximately 
unity at low currents to approximately two at high currents (and 
larger values when base push-out occurs). The shift to a value of 
two represents high injection effects. The n = 1 and n = 2 asymp- 
totes intersect approximately at I c = I k and V eb = V k (see Fig. 2). 



If wc define an emitter capacitance C e by 

(39) 



' dV. b 



then the effective emission coefficient at low current values, where 
the current contribution to Q b , equation (14), is negligible, is 
given by* 

. - 1 + ®^a. ( 40) 

Thus, small deviations from the ideal exponential law are caused 
by the emitter capacitance. These deviations are present even at 



* If the transistor is used in a common emitter configuration it may be useful to 
define the effective emission coefficient as in equation (39), but with Vet instead 
of Vcb held constant. For this case, shown in Fig. 2, the emitter capacitance C« in 
equation (40) is to be replaced by the sum of emitter and collector capacitances. 
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Fig. 2 — Collector and base currents for three values of collector-emitter voltage 
versus base-emitter voltage. Also shown are "knee" point (Vr, Ik), and slopes 
corresponding to values of 1 and 2 for the emission coefficient n. 



low forward currents 1 . Additional increases of n at intermediate 
currents are caused by base resistance (not considered in this 
paper), 
(u) For currents low enough such that base widening effects are 
negligible, the emitter-collector delay time r d for common emitter 
operation is given by 



Td = 



dQj, 





, kT{C. 


+ C e ) 


' ! 


1 q 


I c 


1 


kT{C. 


+ C.) 



(41) 



Q b 



t In principle, equation (40) could be used to obtain the emitter capacitance from 
a dc semilog plot of /,- vs. V.b. However, very accurate temperature control would 
be required. 
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At low current values, and for C e <<C C. , the denominator of equa- 
tion (41) is approximately 1/n, equation (40). Then the value of 
the emitter capacitance C. may be obtained from the slope 
n(kT/q)C. of a plot of delay time r d versus reciprocal collector 
current (see Fig. 6). It is this forward-bias emitter capacitance 
that may be used to set the parameter a, 2 (see equation 1 of Ap- 
pendix C). 

As a demonstration of the proposed model, various characteristics 
for a double diffused silicon npn transistor at room temperature (kT/q = 
0.02585 V) were computed with the parameter values listed in Table 
II. The results are shown in Figs. 2 through 7. We treat the intrinsic 
transistor model, without inclusion of parasitics, such as base resistance 
and inactive base region. Figure 2 is a semilog plot of collector and base 
current vs. emitter-base current for V e , = 1, 2 and 3 V. Note the transi- 
tion from n = 1 to n — 2 near I k = 19 mA and V k = 0.74 volts. Figure 
3 shows common-emitter low-frequency current gain vs. collector 
current for various collector voltages. The frequency dependence is 
conveniently characterized by the low-frequency approximation to the 
unity-gain frequency, 24 j L ■ (For high-current-gain transistors in the 
active region, j L is synonymous with the conventional cut-off frequency 
j T .) Figure 4 shows j L vs. collector current for various collector voltages. 

Table II — Model Parameter Values 

h = -1.875 X 10- 2 amps. 
v k = 28.7 

t, = 3.2 X 10- 10 secs. 
r, = 10.0 

ii - 2.35 X 10- 4 
i 2 = 2.19 X 10-' 
n. = 1.5 
it = 7.65 X 10- 3 
n c = 1.5 

Vo, = 27.1 

m, = 0.24 

a el = 0.337 

ffl.a = 1.03 X 10" 2 

V 0C = 27.1 

m c = 0.1265 

a e i = 0.187 

c c2 = 7.17 X 10-' 

Vrp = 18.2 
r„ = 10.0 
r p = 4.55 
n„ = 3.0 
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Fig. 3 — Small-signal low-frequency current gain h/, 
three values of collector-emitter voltage 



versus collector current for 



Figure 5 presents a family of I c vs. V ee characteristics, with I b as param- 
eter. Figure 6 gives results of simulated r d vs. l/I c measurements of the 
model. The same information is displayed as j L contour plots in Fig. 7. 



v. CONCLUSION 

With the aid of a new charge control relation, a compact bipolar 
transistor model, suitable for network analysis computer programs, has 
been developed. This model is an "integral model" in that parameters 
important for the ac response also shape the dc characteristics. The 
model includes high-injection effects; for low bias and idealized emission 
coefficients it is equivalent to the Ebers-Moll model. Figures 2-7 exhibit 
many features of real transistors that conventional models do not con- 
tain, or can give only by curve-fitting with the use of many parameters. 
These figures were obtained for an implementation of the model using 
21 parameters (Table I, Appendix B). For several of these parameters, 
values can be estimated a priori. Such values may be used as default 
values. Sometimes users may be willing to make some sacrifice in the 
accuracy of the model in return for having to specify only few key 
parameters. It appears feasible to require that only I k (proportional to 
emitter area), the Early voltage V A , maximum /3 (at some collector 
voltage, for example, V ee = 5 V), maximum j L (at V ce = 5 V) and the 
collector current at which maximum f L occurs, need be specified and 
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Fig. 4 — Low-frequency approximation to unity-gain frequency //, versus collector 
current, for three values of collector-emitter voltage. 
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Fig. 5 — Collector current versus collector-emitter voltage for various values of 
base current. 
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Fig. 6 — Emitter-collector delay time versus reciprocal collector current for three 
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0.04 




2 3 4 5 

COLLECTOR EMITTER VOLTAGE V ce IN VOLTS 



Fig. 7 — Contours of constant f L in collector-emitter voltage, collector-current 
plane. 
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that the full set of model parameters is automatically computed with 
the use of default values. 

Not included in this paper, but to be included in careful modeling, 
are parasitics such as base resistance, inactive base region, collector 
resistance, and so on. The base resistance may be modeled as consisting 
of a constant, exterior resistance incurred in the inactive base and an 
interior resistance. To reflect conductivity modulation of the base, 
the value of the interior resistance may be made inversely proportional 
to the base charge Q b . It appears feasible to extend the intrinsic model, 
for example to include carrier generation by impact ionization. 
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APPENDIX A 

Base Push-Out Effect 

The approach towards modeling the base push-out effect was guided 
by results obtained in a detailed analysis of this effect. 25 The discussion 
is based on pnp polarity (I e is a negative quantity). Assuming constant 
resistivity p in the collector region adjacent to the base (epitaxial region), 
the base push-out effect starts approximately at a collector current 
value 

where A, is the emitter area and W c is the width of the lightly doped 
collector region. Let W„ f i be the effective width of the base. For —I e < 
I x , the effective base width is equal to the metallurgical base width, W b , 

TF.ii = W b . (43) 

For — I c > I \ , the effective base width is approximately given by 



-8- 



W mtt = W b + W e \l -^J- (44) 

Equations (43) and (44) can be written as 

WM . Wt + I£UL±l£LzJL±U. m 

Equation (45) gives much too abrupt a transition for W efl from W b to 
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(W b + W e ), and becomes numerically unstable when I e approaches 
zero. Hence, we modify equation (45) with an additional parameter 7 2 , 
such that 

Wm . Wb + k. w. + i,r+ i f-v. + w . (4C) 

The low-current forward transit time 77 is to be modified by the square 
of the ratio of the effective base width to metallurgical base width, to 
give the total base transit time r fb : 

'» = -fe) 5 s " B (47) 

with 

B ~ v + 4Tf t /: + n ( (4bj 

The quantity B may be expressed in terms of model parameters and the 
normalized collector current i e = I c /Ik 

B = ll+^ i^ + 2 +^-(v + ^)l 2 \ 2 (49) 

14 il + r p J 



with 



and 



Wr. I 



Tw w>' Vp \Ik 



. = (V oe - V rH )A e _ v oc - v, - v r 



pWJ k v rp 



(50) 



So far, equation (49) models effects in a one-dimensional transistor. 
Emitter crowding and carrier storage in the inactive base cause the 
transit time at high currents to increase more strongly than given by 
equation (49). For a first-order modeling of emitter crowding we replace 
the exponent 2 outside the square brackets by an adjustable model 
parameter n p ("push-out exponent"): 

The quantities r a , r p , v TP , n p and v oc are model parameters (groups 5 
and 4). 
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APPENDIX It 



A Summary oj the Model 

Terminal quantities are: collector current I c , base current I b , emitter- 
base voltage V cb , collector base voltage V cb , and base charge Q b . 
Define the following normalized quantities: i e = I c /I k , i b = I b /I k , 
v. = (qV eb /kT) - Vk ,v c = (qV cb /kT) - v k , q b - Q b /(I k r t ), and define 
the function 

Then the following relations hold between the normalized quantities 

ib, = iiifi" — e k ) 4- i,(e" /n ° — e ke ), (52) 

ibc = i 3 (e" /n ' — e kc ), (53) 

h, = —ib. - ibc , (54) 

it = i c 4- (v„r - v e — v k )/v rv , (55) 

B -\ l + r » 4 (*1 + r p ) /' ( ' ,(,) 

?, = ! + /(». 4- v k , P.) + f(p c + v k , P.), (57) 

g 2 = B(e v ' - e k ) 4- r,(e° - e k ), (58) 

7. = qJ2 4- [(g,/2) 8 + g 2 ]*, (59) 

*■ = - (C " 'n ^ + *• ■ <«» 

APPENDIX C 

Discussion oj Model Parameters 

The knee parameters 7 t and » 4 , and 77 are related to the conventional 
intercept current I, (saturation current), the base charge at zero bias 
Q bo , and the nominal value of the low frequency approximation to the 
unity-current-gain frequency j L , by 

h = =&* , (61) 



»' = '" (*) • 



(02) 
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and 



Tf = 



2tt/ l 
In terms of structural parameters r, is given approximately by 



(63) 



t/ - srfc (64) 



2 V D 

where w b is the base width and where 77 represents the drift effect in the 
base. 28 ?j is unity for uniform base doping and has typical values between 
two and ten for diffused-base transistors. 

The zero-bias base charge is approximately given by 

Q b0 = -A.qN b (65) 

where A, is the emitter area, q the electronic charge and N b the number 
of impurities per unit area in the base. Typical values for N b are a few 
times 10 12 per cm 2 (lower values would cause premature punch through 
and high values cause low injection efficiency and/or low j L ). The 
intercept current is given by 

N„ Q bo 

where n ( is the intrinsic carrier concentration and D is the effective 
diffusivity of carriers in the base. From equations (62) and (64) through 
(66), one obtains 

2*1^^ + ^277. (67) 

n,- 

The parameters *, , i 2 , H determine the current gain of the transistor. 
The variation of forward gain with current at low currents is governed 
by n, and the inverse gain variation by n c . 

The group 3 parameters model the emitter junction capacitance. The 
offset voltage V„ is approximately the conventional "built-in" voltage, 
which has a typical value near 0.7 volts for silicon at room temperature, 
or v ot = V,./(kT/q) = 27. The grading coefficient m, depends on the 
type of doping transition: it is one-fourth and one-sixth for ideal step 
and linearly graded junctions respectively. Typical values for emitter 
junctions are in the neighborhood of 0.2. Parameter a tl is related to 
the zero-bias capacitance C. by 23 



C .V ot 
Q b .(l - 2m.) 



o.» - „ „~ n_x - (68) 
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If one assumes the emitter junction to be a step junction and the base 
doping to be uniform, and if one expresses base doping and base width 
in terms of the number of impurities per cm 2 in the base, N b , and in 
terms of j L , then an approximate formula for a ei may be derived to be 



(/,/10 9 Hz)* 
" A (N b /5 X 10 12 cm" 2 )* ' UUJ 

where e is the permitivity and D the diffusivity of electrons (holes) 
in an npn (pnp) transistor. For a silicon npn transistor, the numerical 
value for A based on equation (69) is 0.147. For an actual double dif- 
fused transistor of j L = 400 MHz, the value of A obtained from param- 
eter fitting was found to be 0.202. 

The last parameter a e2 in this group is related to the forward bias 
capacitance, C. f , deduced from the slope of delay time versus reciprocal 
emitter current by 23 



a,, = 



pi - ye,] - (71) 



where r is a numerical coefficient approximately equal to unity, the 
exact value depending on the doping profile. Typical values for a, 2 
range from 10 -2 to 10~ 3 . If emitter capacitance effects are not of im- 
portance to the users, and if the user does not want to specify group 3 
parameter values, the following default values are suggested: 

v oe = 27, silicon; 

m t = 0.2; 

(/,/10 9 Hz) 1 



a., = 0.2 



(N h /o X 10 12 cm" 2 ) 5 ' 
a. 2 = 3 X HT 3 . 

The group 4 parameters model the collector junction capacitance. 
The parameters v oe , m e , a el and a e2 have similar meanings as their 
counterparts in the emitter junction capacitance. Typical default values 
(for silicon transistors) for v oe , m e , and a e2 are 

v oc = 27, silicon; 

m c = 0.15; 

a r o = 10" 3 . 
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We shall relate a el to the output characteristics of the transistor. The 
Early voltage as defined in Section II is of magnitude comparable to 
that of the punch-through voltage V T , defined as that voltage for which 
the charge associated with collector capacitance, Q ce , equals minus 
Q bo . We denote the coefficient relating the Early voltage and the punch- 
through voltage by r A 

V T = r A V A . (72) 

Then a e , is given in terms of the Early voltage by 

«.. = T-n ^TTT^ (73) 



rl / ir ir \i-2»ir 

r A V A — V nr \ _ 

\ V oc ) 



1 



The exact value of r A depends on details of the doping profile and on 
the region in the I e vs V ee domain from which the Early voltage is ex- 
trapolated; a typical value of r A is 1.7. Thus, given the Early voltage 
and the other group 4 parameters, equation (73) is convenient for 
estimating a cl . It should be noted that a cl refers to the collector capaci- 
tance of the intrinsic transistor. The terminal collector capacitance 
will be dominated by that of the inactive base region. 

The group 5 parameters model the base push-out effects. V TP = 
(kT/q)v rv is the resistive voltage drop across the collector, caused by a 
current of magnitude I K • The ratio of the width of the collector epitaxial 
region to the width of the metallurgical base is designated r„ . The 
parameter r„ determines the steepness of the variation of the forward 
delay time as a function of a collector current in the current range 
where base push-out is incipient. The base push-out exponent n v deter- 
mines the fall-off of j L for high currents. For n v = 2, j L has a tendency 
to level off after it has decreased from its maximum value by a factor 
of (1 + r, r ) 2 . For n v > 2, the decrease continues beyond this level. In 
the example given in Section IV of this paper, the following values for 
v rn , r w , r p and n p are used. 

v rv =18.2, 

r„ = 10.0, 

r„ = 4.55, 

71 „ = 3.0. 
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